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ABSTRACT 

We calculate the contribution to the cosmic infrared background from very massive 
metal- free stars at high redshift. We explore two plausible star-formation models and 
two limiting cases for the reprocessing of the ionizing stellar emission. We find that 
Population III stars may contribute significantly to the cosmic near-infrared back- 
ground if the following conditions are met: (i) The first stars were massive, with 
M > 100 Mq. (ii) Molecular hydrogen can cool baryons in low-mass haloes, (iii) Pop 
III star formation is ongoing, and not shut off through negative feedback effects, (iv) 
Virialized haloes form stars at about 40 per cent efficiency up to the redshift of reion- 
ization, z 7. (v) The escape fraction of the ionizing radiation into the intergalactic 
medium is small, (vi) Nearly all of the stars end up in massive black holes without 
contributing to the metal enrichment of the Universe. 

Key words: cosmology: theory - early Universe - galaxies: formation - stars: for- 
mation - diffuse radiation - infrared: general. 



1 INTRODUCTION 



It has long been realized that observations of the cosmic in- 
frared background (CIRB) can place important constraints 
on the energy production associated with the formation 
of cosmological structure (see Hauser & Dwek 2001 and 
references therein). The cosmic background is the locally 
measured radiation density from all extragalactic sources. 
Starlight dominates the CIRB in the near-IR, whereas the 
mid- and far-IR backgrounds result primarily from dust 
emission (e.g., Dwek et al. 1998). 

Observational efforts to measure the near-IR CIRB are 
hampered by strong and uncertain foreground contamina- 
tion from zodiacal dust. Cambresy et al. (2001), using data 
from the Diffuse Infrared Background Experiment (DIRBE) 
on the COBE satellite, found that the integrated light from 
all galaxies (estimated from deep surveys) is not sufficient 
to account for the measured CIRB at 1.25 /.im (J band) or 
2.2 fim {K band). Wright & Johnson (2001) analyzed the 
same DIRBE data as Cambresy et al. (2001), but subtracted 
a different zodiacal light model; they find a J-band CIRB 
that is consistent with the contribution from galaxies, but 
their JsT-band CIRB value is larger than the background in- 
ferred from galaxy counts. Totani et al. (2001) model the 
contribution from galaxies missed by deep galaxy surveys 
and find that it is unlikely to be greater than 30 per cent of 
the measured light. 

Bond, Carr & Arnett (1986) suggested that the first 



generation of stars in the Universe, called Population III 
(Pop III) stars because they are assumed to have near-zero 
metallicity, may contribute significantly to the cosmic back- 
ground. The cosmological impact of Pop III stellar radia- 
tion has been addressed before (e.g., Carr, Bond & Arnett 
1984; Bond et al. 1986; Couchman & Rees 1986; Haiman & 
Loeb 1997; Ciardi et al. 2001). Recent theoretical advances 
have improved our understanding of the physics of Pop III 
star formation (see Barkana & Loeb 2001 and references 
therein); in particular, three-dimensional numerical simula- 
tions suggest that the primordial initial mass function (IMF) 
may have been concentrated at stellar masses ^ 100 Mq 
(Bromm, Coppi & Larson 1999, 2002; Abel, Bryan & Nor- 
man 2000). Sta rs with masses ^ 100 Mq, which we will call 
'very massive' (Carr et al. 1984), have spectra and luminosi- 
ties that asymptotically approach the blackbody form and 
the Eddington limit, respectively (e.g., Bromm, Kudritzki & 
Loeb 2001b). 

In this paper, we investigate whether a significant frac- 
tion of the near-IR CIRB may come from very massive 
Pop III stars. Since these stars are very luminous per unit 
stellar mass over a short lifetime, they may contribute sub- 
stantially to the CIRB. We use two different models for the 
formation of dark matter (DM) haloes to calculate the rate 
at which baryons are processed through haloes with virial 
temperatures high enough to allow baryonic cooling {^). In 
addition, we consider two possibilities for the re-processing 
of the stellar emission by the gas in the halo and by the 
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intergalactic medium (§^). In §^we show our model results 
and the observational data. Section ^ contains a discussion 
of our results and compares them with the data. 

We adopt the following values for the cosmological pa- 
rameters: h = 0.7, D.m = 0.3, = 0.7, Qsh'^ = 0.019, 
(78 = 0.9. Here h is the dimensionless Hubble constant, 
-ffo ~ lOO/i km s~^Mpc~^; Q,m, and JIb are the total 
matter, cosmological constant and baryon density in units 
of the critical density and as gives the normalization of the 
power spectrum on the 8ft~^Mpc scale. 

2 STAR-FORMATION RATE 

2.1 Overview of the star- formation model 

We assume that the star-formation rate (SFR) of Pop HI 
stars is related to the cooling of baryons in collapsed dark- 
matter haloes (henceforth just 'haloes'). The cooling rate 
of baryons is a very strong function of baryon temperature 
near certain temperature thresholds: above these thresholds 
a new cooling mechanism can act within the gas. For a given 
choice of cooling mechanism, we would like to calculate the 
rate at which baryons are heated above the associated tem- 
perature threshold, and thus may collapse to high density 
and form stars. 

If baryons in a halo are shock-heated to the virial tem- 
perature of their halo, the temperature of those baryons is 
related to the mass of that halo according to T cx M^^^ . So 
as a proxy for the temperature history of the baryons, we 
may calculate the mass history of haloes, assuming a con- 
stant ratio of baryonic matter to dark matter in all haloes, 
and that baryons below the cooling threshold are always 
heated to the virial temperature of their halo. In this pic- 
ture haloes above a critical mass may form stars, because 
those haloes have heated their gas above the critical tem- 
perature, and haloes below the critical mass do not form 
any stars. Additionally, we parametrize the fraction of eligi- 
ble baryons that actually do form stars with a constant star 
formation efficiency, 77. 

The mass-assembly history of dark-matter haloes is 
computed using the extended Press-Schechter formalism 
(Lacey & Cole 1993), where the a{M)-M relation is eval- 
uated with the power spectrum of Eisenstein & Hu (1999). 
Here a{M) is the standard deviation of the linear density 
field smoothed on scales containing a m ean mass AI. Our 
choices of Merit (z) are described in §2.3, 



We consider two models for the SFR. In both mod- 
els, star formation is triggered in a halo as it accumulates 
enough mass to put it above the threshold mass Merit (2:), 
but the models differ in how they treat additional material 
that merges into such haloes. In the 'ongoing' model, all gas 
in haloes with M > Merit (2) is eligible for star formation. 
Conversely, the 'single-burst' model only allows star forma- 
tion in haloes that have not previously formed stars and have 
no progenitor halo that has formed stars, i.e., both merging 
haloes are crossing the critical threshold for the first time. 

2.2 Specific star-formation models 

2.2.1 Ongoing star formation 

In the ongoing model, the SFR in a single halo is propor- 
tional to the growth (in mass) of that halo, after it has be- 



come more massive than Merit (z). Star formation does not 
occur in a quiesciently evolving (i.e., non-accreting) halo; 
the star formation is not continuous in the usual sense. But 
star formation is not inhibited after the first generation of 
Pop III stars is born in a halo, which assumes that radi- 
ation and mechanical outflows from the life and death of 
the star(s) have no impact (zero feedback) on future star 
formation. 

The SFR in the ongoing model is proportional to the 
time derivative of the total mass contained in haloes that 
have mass M > Merit (z). The SFR per comoving volume, 

V'on, is 



i^oniz, Mcvlt{z)) = T] 



Ob d 

fim dt 



dMM^(M,.). (1) 



Here ups = nps(M, 2) is the comoving number density of 
haloes of mass M at redshift z, given by Press & Schechter 
(1974). The integral expresses the collapsed mass per co- 



moving volume contained in haloes above the critical mass. 
The time derivative converts this to a mass rate, and the 
prefactors convert from total mass to stellar mass. 



2.2.2 Single-burst star formation 

In the single-burst model, a halo forms stars when the ac- 
cretion of matter pushes its mass above the critical mass. 
Merit (2); that is the only time that the halo will form stars, 
regardless of additional mergers and mass accretion. This is 
a very extreme model which assumes Pop III stars make a 
long-term change in their surroundings which permanently 
prohibits future Pop III stars from forming there. This ex- 
treme model could be realistic if Pop HI stars heavily enrich 
the interstellar medium with metals. Then, though stars will 
likely eventually form out of that material, Pop III stars 
will not. The SFR may lie somewhere between the pre- 
diction of this model, and the prediction of the ongoing 
model discussed above (in which all baryons in haloes with 
M > Merit (z) are eligible to form stars). 

The SFR in the single-burst model can be calculated 
with the extended Press-Schechter formalism: A halo grows 
in mass through discrete mergers and gradual accretion, 
which is treated as a series of small mergers in extended 
Press-Schechter theory. As an example, take a halo with 
mass Ml at some redshift 2 that merges with mass AM for 
a total mass of M2 — Mi + AM at redshift 2+d2, where d2 is 
small and negative, i.e., a short time later. If Ali < Merit (z) 
and il/2 > Merit (z), then, provided AM < Merit (2), the 
halo undergoes a burst of star formation. We parametrize 
the mass of stars formed, M*, by 



M,=n^M2. 

The SFR per comoving volume, Vburst, is 



(2) 



V'burst(z, Merit (2)) 



dMi 



dni 



dMi 



Af„it(z) + Mi 



dM2 M2 



1 JM^itM 



AM2 dt' 



(3) 



where ni = nps(Mi,2) and P = P{Mi, M2, z) is the prob- 
ability that a halo with mass Mi merges to a new mass 
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M2 > Ml at redshift z ( |Lacey fc Cole 1993| ). The pref- 
actor of 1/2 corrects for double counting in the integrals 
(for fixed Mi and AM, the integral as written counts both 
Mi + AM = M2 and AM + Mi = M2).0A similar approach 
to model the history of star formation at high redshifts has 
been discussed by Barkana & Loeb (2000). 



2.3 Critical halu maws for star furiiiaLiuii 

There is a relationship between the mass of a halo and its 
virial temperature. Using the assumptions we made in § |2.l[ 
we can substitute a critical temperature for baryonic cooling, 
Tcrit, into that relation to find the corresponding critical 
mass: 



Merit (.z) = 0.94 X 10** Mq 



f — Y ( — 
V0.7/ VO.J 



-1/2 



V 10 y Vo.ey uo* k) 



(4) 



where fi is the mean molecular weight (/i = 0.6 for ionized 
gas [Tcri t <: 10"^ K1 and m =, 1-2 for neutral gas [Tcrit ^ 



10" K]) (Barkana fc Loeb 2001) 



We seek to model only two cooling processes: radiative 
cooling from molecular hydrogen and radiative cooling from 
atomic hydrogen. For molecular cooling, we choose Tcrit = 
400 K (e.g., Tegmark et al. 1997; Abel et al. 1998). For 
atom ic cooling, we choose Tcrit ~ 10* K (Barkana & Loeb 
200l|). 



2.4 Limits on tiie abundance of Pop III stars 

A realistic model of Pop III star formation must conform 
to some observational constraints. The matter density pro- 
cessed through Pop III stars as a fraction of the critical den- 
sity is 



f^iii(z) = — 

pc 



dt' 



dz' 



dz' , 



(5) 



= [Ho{l + z)E{z)]- 



with 
dt_ 
di 

and, in a flat universe, 

E{z) = [Q^{l + zf + Qa]'^^. 

Here, ip is the SFR and pc is the density required to close the 
Universe. A strong upper limit on Qm is that it must not 
exceed Ob . Below we discuss two constraints on the possible 
products of Pop III stars, black holes and heavy elements 
(see also Schneider et al. 2002). 



2.4-1 Black holes 

A nonrotating star with mass above ~ 260 Mq is predicted 
to ev olve directly to a black hole w ithout any metal ejec- 
tion (Fryer, Woosley & Heger 2001). The fate of rotating 



stars in that mass range is unknown, but may be similar. 



* There is some ambiguity in the SFR arising from the extended 
Press-Schechter merger rates. We defer discussion of that to future 
work, and here evaluate the SFR as written. 



If so, and if the Pop III IMF were concentrated at masses 
above ~ 260 Mq, most of the baryons once in Pop III stars 
would now be in black holes. Such black holes would popu- 
late the halo of our galaxy, but they would be very difficult 
to detect unless they have a companion; e.g., the MACHO 
microlensing study places no constraint on the fraction of 
the Milky Way halo contained in ^ 100 M© black holes (IaT 



cock et al. 2001). Detailed theoretical studies may eventually 



reveal the fraction of such black holes that currently have 
a companion object, and searches for those binary systems, 
or future microlensing studies, could lead to constraints on 
SIbh (e.g., Agol & Kamionkowski 2001). But currently even 
if f^BH = fiiii there is no additional constraint on Qui from 
limits on black holes. 



2.4-2 Metal enrichment 

A nonrotating star with mass between ~ 140 Mq and ^ 
260 Mq is expected to end its evolution in a pair-instability 
supernova , which completely disrup ts the star and leaves no 
remnant (Heger & Woosley 2002). In this case all of the 



nucleosynthetic products of the star are ejected, and poten- 
tially pollute the intergalactic medium (IGM) (cf. §2.4.1 



where all were swallowed into a black hole). Thus metal 
abundances in the Universe, especially at high redshift, are 
limits on the number of Pop HI stars which resulted in su- 
pernovae, given assumptions about the mixing of the metals. 

If a fixed fraction tpi ('pi' for 'pair instability') of 
the mass converted into Pop HI stars is expelled in pair- 
instability supernovae, the nucleosynthetic yield calculation 
from Heger & Woosley (2002) may be used to determine 
the contribution from Pop HI stars to the metallicity of the 
IGM. For 



log 



f? 



(6) 



where i is a species and and /® are the mass fractions of 
species i in the Universe and Sun, respectively. Similarly, if 
we denote the ejected mass fraction of species i by , then 



and 



^^III rE 

^piT;-.A 



log 



(.. 



-P 



(7) 



(8) 



where Pi = jf I jf is the production factor of species i. 
Heger & Woosley (2002, table 4) calculated Pi for many 
isotopes and stellar masses. We here take ~ 200 Mq as a 
fiducial mass for a pair-instability supernova progenitor. 

A few values of interest from their table are Pi(i = 
i2(^ i6q 24^g 28gi 56pg) = (13.2,45.8,85.7,353,49.8). In 
general, the highest values of Pi are for the alpha elements. 

Observations of the abundance of species i in the high- 
redshift universe are sometimes quoted in slightly different 
notation from [i]. Assuming epifiin <C Ob, [j/H] — [i]. Us- 
ing tabulated values for (Heger & Woosley 2002, table 
3), it is straightforward to calculate Q.i — J^Qb- With as- 
sumptions about the mixing of the ejecta from Pop HI stars, 
measurements of metals in the IGM at high redshifts (or the 
metallicities of stars formed at high redshifts) can place a 
limit on a combination of the IMF of Pop HI stars (which 
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Figure 1. The star-formation rate (SFR) for tlie ongoing and 
single-burst models, for critical temperatures of Tc = 400 (solid 
lines) and 10* K (dashed lines). For each the upper curve is 
the ongoing SFR and the lower curve is the single-burst SFR. 
All curves are for r] = 0.4; each curve is labelled by the value of 
/b = Oiii(2 = 7)/Qb it implies. 



Figure 2. The matter density processed through Pop III stars 
in units of the critical density, as a function of redshift. The four 
curves correspond to the four different SFRs of Fig. |^. The hori- 
zontal dashed lines labelled Hi, He I and Hen correspond to the 
Qui required for Pop HI stars to produce 10 ionizing photons 
per particle of each species. These lines are only relevant for the 
/esc = 1 case; see section H. 



we simplify into the parameter epi) and Qm. 

-3 



For example, 



if Epi = 1 and /b = ^Iui/Qb — 3 x 10"'', then Pop III stars 
would, in the mean, enrich the Universe to the solar abun- 
dance of silicon already at a redshift ^ 7. In contrast, if 
Epi — 0, then metallicity measurements place no constraints 
on the abundance of Pop III stars. 



2.5 Star- formation model results 

To compare our models, we present the SFR as a function 
of redshift in Fig. |l| We assume that Pop III star formation 
shuts off at 2end — 7; this may be due to photoevaporation o f 
low-mass haloes due to reionization ( Barkana & Loeb 199!: ) 
or some other process. Figure ^ shows the SFRs normalized 
to r; = 0.4. The curves are labelled with the corresponding 
value of /b. 

Figure I shows niii(z) for the SFRs of Figure |l|. The 
horizontal lines are explained in §n. 



3 SPECTRA 

3.1 Input stellar spectrum 

Metal-free very massive stars have spectra simila r to a 
~ 10^ K blackbody spectrum (|Bromm et al. 2001b|). Most 



of the energy is then radiated in photons with energies 
> 13.6 eV. There are also many photons produced capa- 
ble of ionizing Hen (but see Schaerer 2002). The specific 
luminosity per solar mass o f very massive stars is almost 
independent of stellar mass (Bromm et al. 2001b). We will 



therefore take the spectrum of a 1000 Mq star as our fiducial 
input Pop III stellar spectrum. 

When a Pop III star becomes luminous, we expect the 
host halo to contain gas not incorporated into stars; we call 
this gas the 'nebula.' We refer to the gas outside the col- 
lapsed halo as the IGM. The physical environment of Pop III 
star formation will likely undergo a transition from dense gas 
in the nebula to a diffuse IGM, but we will model all gas as 
either part of the nebula or part of the IGM. We also assume 
that both the nebula and IGM are entirely free of dust. 

Since the nebula and IGM are neutral in the absence 
of stellar radiation, both have the potential to play an im- 
portant role in reprocessing ionizing photons from a star. 
We examine two limiting cases for the importance of each of 
these phases: In the first case, Pop III stars are enshrouded 
in dense nebulae, and all of the reprocessing of ionizing ra- 
diation takes place in the halo (the IGM still plays a role in 
scattering Lya photons). In this case, the escape fraction of 
ionizing photons from the nebula, fesc, is zero. In the second 
case, the nebula plays no role and fesc ~ 1; all reprocessing 
occurs in the IGM. In the rest of this section, we discuss the 
reprocessed spectrum of Pop III stellar radiation for each of 
these cases. 



3.2 No escape of ionizing radiation into the IGM 

3.2.1 Properties of a nebula 

Numerical simulations suggest that when a Pop III star 
forms, the nebula consists of a higher density phase, with 
riH — 10* cm"'', and a lower density phase (e.g., Bromm et 
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al. 1999). We make the simplifying assumption that half of 
the nebula's mass is contained in a homogeneous phase with 
density nn = 10* cm~^ that completely covers the star(s), 
and ignore the lower density gas. We take the mass frac- 
tions of hydrogen and helium as X = 0.75 and Y = 0.25, 
respectively. 

Ionizing radiation from the star(s) creates an H ll re- 
gion in the dense nebula. Because of the hardness of our 
input spectrum, in the inner part of the Hll region helium 
is doubly ionized (the Hem region). In the outer part it is 
singly ionized (the He ll region; the spectrum is hard enough 
that there is no Hll/Hel region). The majority of photons 
above the He ll ionization threshold ionize He ll rather than 
Hi. We iteratively solve for the sizes and temperatures of 
these regions, using the thermodynamic equations from Cen 
(1992). The Hem region, comprising 0.4 of the volume of 
the Hll region, has a temperature of 3.6 x 10* K; the Hell 
region is cooler, at 2.7 x 10* K. In both regions the pri- 
mary cooling mechanism is Hi recombination, but free-free 
emission and cooling via collisional excitation of H I are also 
important. The total nebular emission is not very sensitive 
to the relative sizes of the Hell and Hem regions. 

Given these properties of an ionized region, we can de- 
termine /osc.^ The volume and mass of the Hem region are 

QHcII 



VhoIII = 



neriHemaB (He ll, T = 3.6 x 10^ K) ' 



Mneiii = A"^H (riH + nne) VhcIii — 0.7 



(10) 



where QhcII is the stellar emission rate of photons energetic 
enough to ionize He ll, QB(He ll, T) is the Case B recombina- 
tion coefficient for Hell, here M» is the mass of the ionizing 
Pop HI star(s) in the nebula and ^ = 1.2 is the mean molec- 
ular weight. Because in the Hem region recombinations to 
Hell provide enough photons to keep the hydrogen ionized, 
the volume and mass of the He ll region are 



BnmiaB (H I, T = 2.7 x lO-* K) ' 



MhcII = M"^H {n-a + uho) VhcII 1.0 ( — ^ 

1V10 



(11) 



(12) 



where Qhi is the rate of H l-ionizing photons, and ae (H I, T) 
is the hydrogen recombination coefficient. 

The total mass of ionized gas is then 1.7 times the mass 
of the ionizing star(s). In the single-burst star-formation 
model, the nebula mass in tih ~ 10* cm~'^ gas, A/ncb, is 



1 1 

2 ~ 



V 



M,. 



(13) 



For 77 = 0.4, Mncb = 0.75 M,, which implies /esc = would 
be difficult to achieve for our model nebula. In the ongoing 
model, though, Mncb > (1 — ri)M,/{2rj) in general, and es- 
pecially at lower redshift. This is due to the short lifetime 
of Pop HI stars: star formation in a halo subsequent to the 
first episode usually takes place after previous generations 



t The dynamical evolution of the gas on the timescale of the 
lifetime of a Pop III star may be important, as the pressure of 
the nebular gas may exceed the gravitational pressure, but here 
we treat the nebula as static. 



71 



O 




Figure 3. Spectrum of a 2 = 15 Po p III star plus nebular emis- 
sion, for the /esc = case of section 3.2, The long-dashed line is 



the spectrum of the star, cut off for absorption shortward of the 
Lyman limit. The dotted line is the spectrum of free-free emis- 
sion from the nebula. The short-dashed line is emission from hya 
recombination in the nebula at z = 15, corrected for scattering in 
the IGM. The solid line is the sum of the spectra. All spectra are 
in the rest frame of the star. 



of stars in the halo have stopped radiating. Thus /c. 
may be possible in these haloes. 



3.2.2 Resulting spectrum 

3.2.2.1 Recombination emission. Approximately 1/2 
of the energy radiated by a Pop HI star in a nebula is ul- 
timately reradiated by recombinations. The mean energy of 
a free electron just before it recombines with a proton may 
be estimated from 



/3B(Hi,r) kT 
aB(Hl,T) ' 



(14) 



where /3B(Hl,r) is the recombination emissio n coefficient 



Using T ^ 3 X 10* K, we find (E) = 1.4 eV ( ^eaton 195^ 



Cen 1992). The average energy radiated per recombination, 
then, is 15.0 eV. Under the on-the-spot approximation, all 
Lyman series photons degrade into n = 2 — > 1 transitions, 
where n is the energy level of an excited hydrogen atom. 
Since 10.2 eV is released per n = 2 ^ 1 transition, about 
2/3 of the recombination energy is released in Lya or two- 
photon emission, and 1/3 of the recombination energy (and 
about 1/6 of the total emitted energy) is emitted in other 
recombination lines and free-bound continuous emission. We 



ignore those processes here; see Schaerer ( 2002 ) for a more 



complete computation of nebular emission lines and free- 
bound spectra. 

The relative importance of Lya emission compared to 
two-photon emission is determined by the effective recom- 
bination coefficients to the 2p and 2s states, respectively, 
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and the coUisional excitation rate from the 2s state to the 
2p state. At T = 3 X 10'* K, 0. 75 Lya photons are emit- 
ted for every H I recombination (Osterbrock 1989; Storey & 
Hummer 1995); we ignore the contribution of two-photon 



emission (about 12 per cent of the total emission) to the 
spectrum. 

Since Lya photons resonantly scatter in neutral hydro- 
gen, they will not travel far in the IGM until their frequen- 
cies are shifted away from the resonant frequency, z^Lya . Pho- 
tons initially scattered blueward of the line resonance will 
eventually cosmologically redshift back into the resonance. 
The result is that the IGM ultimately scatters all Lya pho- 
tons to the red side of the line resonance (broadening from 
motions inside a halo doesn't contribute to the line profile). 
For a homogeneous, expanding IGM, the resulting sc attere d 
line profile, (^(z/, z), was simulated by Loeb & Rybicki (1999), 
and we fit their result with 



ly, (z) y-^ exp [^^^^] if > 



if < 



v.^{z) = 1.5 X 10" Hz 



0.019 



(^) 



-'{l + zf 
E{z) 



(15) 



(16) 



This profile results in a strong, asymmetric Lya emission 
line near 1220 — 1225 A with a scattering tail extending to 
long wavelengths. 

In the He ill region. He ll recombinations sometimes pro- 
duce more than one photon capable of ionizing Hi. Recom- 
binations directly to the n = 2 state produce Hell Balmer 
continuum photons, which are capable of ionizing Hi. In 
addition, two-photon decay fro m the 2s state pr oduces 1.42 
H l-ionizing photons per decay ( Osterbrock 1989 ) . Since He ll 
Lya is also capable of ionizing hydrogen, the mean number 
of hydrogen ionizations per He ll recombination is 



= 1.7, (17) 



02 (He II, T) + 1.42alf (Heii,r) +al^(Heii,r) 
QB(Hell,r) 

where 02 (He 11, T) is the recombination coefficient for re- 
combinations directly to the n = 2 state and ot^ (He 11, T) 
i s the effective total recom bination coefficient to the 21 state 
( ^torey fc Hummer 199^ ). 



3.2.2.2 Free-free emission. Free-free radiation ac- 
counts for about 1/4 of the cooling in the nebula. This en- 
ergy is radiated in a continuous spectrum. 



jf = 7.2 X 10-'" "'^'^^ (I) 
z 



-1/2 



X exp 



ergs ^cm^ Hz ^ster ^, (18) 



where is the specific emission coefficient, v is the fre- 
quency of emitted radiati on and nz is the number density 
of ions of net charge Z (Ferland 198C) A Gaunt factor of 
1.3 has been assumed; this results in an error of less than 8 



per c ent over optical and UV frequencies (Karzas & Latter 



1961) 



The luminosity per solar mass of Pop HI stars from free- 
free radiation, is 



iff 47r / .fl HcIII , .fr.Hell \ 



(19) 



Here He ill and He 11 label the emission coefficients and vol- 
umes computed for the Hem and Hell regions; in each case 
the free-free emission is dominated by Hll. 

3.2.2.3 Emission from collisional excitation of Hi. 

Collisions between a free electron and trace H I atom some- 
times excite the H I atom, which then radiates away the 
excitation energy. The collisional excitation rate coefficient 
for transitions from the n — 1 state to state u, q\u, is 



8.629 X 10"^ Q.il,u) 
qiu = =7rT5 exp 



yl/2 



LUl 



kT 



(20) 



where T is in K, i}{l,u) is the (temperature-dependent) ef- 
fective collision strength for transitions from the n = 1 to 
state u, LJi is the statistical weight of the n = 1 state and 
x{ l,u) is the energy difference between the n — 1 and state 
u dOsterbrock 1989[). 



We compute the collisional excitation to the n = 2 and 
n = 3 states using cross-section s from Callaway ( 1985 ) and 
Callaway, Unnikrishnan & Oza (1987). These excitations re- 



sult in additional Lya emission: excitations to the 2p, 3s and 
3d states radiatively decay to n = 1 via Lya photons, and 
atoms in the 2s state (resulting from collisions to either the 
2s or 3p state) may be additionally collisionally excited to 
the 2p state. The rate of collisional Lya emission per unit 
stellar mass, gLyai is 



coll cff 
Qhya = rtcUm gi,2p 



where 

<li%p = ^ L{u) qiu 

U 

and 



L{u) = 



1 



if « G (2p, 3s, 3d) 



0.33 ff«e(2s,3p) 



(21) 



(22) 



(23) 



The factor of 0.33 acc ounts for 2s to 2p collisional excita- 
tion (Osterbrock 1989). Summing over the Hell and Hem 



regions, = 3.4 x lO^" 



3.2.2.4 Total spectrum. For the /esc = case, the to- 
tal specific luminosity per unit stellar mass emitted from a 
Pop III star and nebula and scattered in the IGM, l1{z), 
has three components that we treat: the truncated stellar 
spectrum, Z*^, the free-free spectrum, and the scattered 
Lya spectrum, l]y'°'{z): 



ll{z)=C + f,+ll^'^{z). 

The truncated stellar spectrum is 



ll if hv < 13.6 eV 
if hi/ > 13.6 eV 



(24) 



(25) 



The Z*^ spectrum is slightly modified between 912 and 
1216 A by scattering in the IGM, when the photons in that 
range are cosmologically redshifted into the Lya resonance 
(e.g., Peebles 1993). This effect is small when 1.75 iy,( z) <g; 



8.22 X 10 Hz, i. e., the width of the scattered Lya line (Loeb 



&Z. Rybicki 1999) is small compared to the frequency differ- 
ence between 912 A and 1216 A. Since 1.75i/*(2 = 30) = 
8.27 X 10" Hz, we ignore this correction to the spectrum. 

© 2001 RAS, MNRAS 000, |-?? 



The Contribution of the First Stars to the Cosmic Infrared Background 7 



T3 



X 

-o 



c 




15 20 25 

Redshift (z) 

Figure 4. Properties of the ionized IGM near a Pop III 
star. There are sets of three curves for five values of Zi = 
(10, 15, 20, 25, 30). In each set of curves, the solid line is the tem- 
perature as a fraction of the initial temperature, the dashed line 
is the ionization fraction and the dotted line is [dx/d^|, followed 
from Zi to z = 7. 



The free-free spectrum is given by eq. (^). The Lya 
spectrum is 

l]^^°'{z) = QLya hULya 0(H.yQ - V,z), (26) 

where gLya is the rate of Lya photons produced per solar 
mass of the ionizing star. That rate is 



qLya = 0.75 {qm — (JhcII + l.Tgncii) + gEy 



(27) 



where qm = Qm/M, and qhcII = QhcIi/M, . The factor 
0.75 represents the fraction of hydrogen recombinations that 
resuh in Lya photons, and the factor of 1.7 accounts for 
the number of hydrogen ionization per Hell recombination, 
computed above. 

Figure I show Z^, lu, and l°{z) for z = 15. In 

this model no Pop III ionizing photons escape to the IGM, 
thus Pop III stars don't contribute to reionization of the 
Universe (cf. section p.3|). 



3.3 Complete escape of ionizing radiation into the 
IGM 

3.3.1 Properties of the IGM 

The second case we consider is that the nebula plays no role 
in reprocessing ionizing radiation from a Pop III star, /esc = 
1. This m ay be because, in contrast to the assumptions we 
made in 



i.2.1, the density of the nebula is low, or because 



the nebula is clumped into high-density regions with a small 
covering fraction, or because the nebula was blown away 
by the star(s). Because the timescales of important IGM 
processes extend beyond the lifetime of a Pop III star, in this 



section we will name the redshift of formation of a Pop III 
star Zi, and then describe the evolution as a function of z. 

We assume the IGM is uniform with baryon density 
niGM(z) = 1.7 X 10~^(1 + z)^ cm~^ (which ignores the 
small fraction of baryons in collapsed haloes), X = 0.75 and 
Y — 0.25. Ionizing photons from a Pop III star stream into 
the IGM and form an ionized region. Because the density is 
low, for Zi ^ 30 recombinations are of little importance on 
the timescale of the star's lifetime, r ~ 2 x lO" yr. For the 
purposes of calculating the properties of the ionized region 
of the IGM, we assume that all of the ionizations occur im- 
mediately; at worst this contributes less than a 3 per cent 
error to the computed spectrum. 

Hem recombines more quickly than Hll, and He I re- 
combinations occur o n a com parable timescale to those of 
Hi. As discussed in j j3.2.2.1 , Hell recombinations produce 
1.7 H I- ionizing photons on average. In the low-d ensity limit, 
a He I recombination ionizes 1 hydrogen atom (Osterbrock 



198!:). We treat He I and Hell ionizations as an extra 1 and 
1.7 Hi ionizations, respectively, and track only hydrogen re- 
combinations. 

For Zi ^ 30, a Pop HI star ionizes a volume of the IGM, 
V(z), initially given by 



V{zi) = 



(gHi + 0.7QHcii)r 2.4x10* 



niGM(z) 



{1 + Zi)3 



The initial temperature in this volume is determined from 
photoionization heating to be Ti ~ 6 x 10* K. Given Zi, we 
solve the temperature evolution of the ionized IGM region 
as a function of z until z = 7, when we assume that the Uni- 
verse reionizes. If Pop III stars reionize the Universe before 
z = 7, their ionized nebulae will have overlapped, and the 
analysis of this section will no longer be appropriate. The 
post-reionization contribution of the IGM to the CIRB de- 
pends on the temperature and ionization state of the IGM, 
but is expected to be small under typical assumptions. But 
the most important effect that reionization would have on 
the Pop HI contribution the CIRB m ay b e the effect reion- 
ization has on the SFR. As noted in j j2.5| , reionization may 
stop or severely curtail the formation of Pop HI stars. 

To compute the temperature of the ionized IGM around 
a Pop HI star as a function of redshift, we set the initial 
temperature to Ti = 6 x 10* K and follow the evolution of 
temperature considering cooling from Compton scattering of 
CMBR photons off of electrons, adiabatic cooling of the IGM 
from the expansion of the Universe, cooling from collisional 
excitation of H I, free- free cooling and recombination cooling 
( |Cen 1992 ). Since all of those processes except adiabatic 



cooling depend on the ionization fraction, x, it is solved for 
simultaneously, using the approximation x ~ Uc/niGM'. 



dx 
dt 



■ QB(Hl,r)niGMa:: 



(29) 



(e.g., Peebles 1993). Figuregplots T(z), x{z) and \Ax{z)ldz\ 
for five values of Zi. 

Initially nearly all the hydrogen is ionized, and Comp- 
ton cooling briefly dominates the thermodynamics. Once a 



small fraction, 



10" 



of the hydrogen atoms have recom- 



bined, cooling via collisional excitation becomes the domi- 
nant process. As the temperature decreases, the recombina- 
tion rate increases, which in turn increases the neutral frac- 
tion and thus the collisional cooling rate. When the tempera- 
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Figure 5. Spectrum of a 2 = 15 Pop II I st ar plus emission from 
the IGM, for the /esc = 1 case of section 3.3 The long-dashed line 



is the spectrum of the star, cut off for absorption shortward of the 
Lyman limit. The dotted line is the spectrum of free-free emis- 
sion from the IGM. The short-dashed line is emission from Ly« 
recombination in the IGM, corrected for scattering. The free-free 
and Lya spectra are integrated from 2 = 15 to 2 = 7, and divided 
by the lifetime of a Pop III star to give a useful normalization. 
The solid line is the sum of the spectra. The sharp peak at 1216 A 
is coUisionally excited hya emission; the broad peak at 1000 A is 
hya from recombinations. All spectra are in the rest frame of the 
star (see Section 3.3.2| for explanation). 



(Ferland 1980). The specific luminosity per unit stellar mass 
is 

y\ -1/2 

X lU" x~ (I + z)" ( 

h V ^ 



ll{z, zi) = 47r j5 = 1.6 X 10'* x^'il + z 



(^) erg s-^ Hz-^Me\ (32) 



X exp 



where the second equality follows from equation (|28[). The 
Zi dependence of results from the implicit dependence of 
T and x on Zi and z. 

3.3.2.3 Emission from collisional excitation of Hi. 

The rate of Lya photons produced by collisional excitations 
is given by eqs. ( pl| ) and (p8[): 

^ /I ™\ I „cff ^™~3 i\/T~i (33) 



= 2.0 X 10^^ :r (1 - a;) (1 -f zf ql%^ cm'^ M"^ 



^0 



In the low density limit, though, collisions from the 2s to 2p 
state are unimportant, so 



1 if M G (2p,3s,3d) 
ifMG(2s,3p) 



(34) 



Additionally, because the temperature in the IGM varies, we 
make the following fits to the data of Callaway et al. (1987): 



n{l, 2p) ~ -2.41 X 10-^ (^) + 0.148 (^) + 0.170(35) 
and 

n(l,3s) + f7(l,3d) 

~ -2.29x10'=^ Vo.0299 (^) +0.116.(36) 

The fits are good to better than 5 per cent for temperatures 
from 1.6 X 10* to 6 x 10* K, where almost all of the collisional 
excitation occurs. 



ture reaches about 1.5 x 10* K, Compton cooling dominates 
again. For lower values of 21, adiabatic cooling eventually 
dominates; for higher values of Zi, recombination cooling be- 
comes important. Free-free emission is important enough to 
be included. 



3.3.2 Resulting spectrum 



3.3.2.1 Recombination emission. As in §3.2.2.1, we 
consider only the Lya component of the recombination spec- 
trum. The rate of production of Lya photons from recom- 
binations, qi^ya, is 



^L^ya = 0.63 (gni - (Jhcii + l.TgHcii) t 



Ax 
lit 



(30) 



where 0.63 is a representative value for the fraction of hy- 
drogen recombinations resulting in Lya emission near the 
peak of the recombination rate. 



3.3.2.2 Free-free emission. The free-free spectrum is 

(I)" 



j« = 7.2 X 10-3^ n?GM ^'^ 



X exp ( ^ y ^ ergs ^cm'^ Hz ^ster ^ (31) 



3.3.2.4 Total spectrum. For the /esc = 1 case, the to- 
tal specific luminosity per unit stellar mass emitted from 
a Pop III star and IGM, including scattering, is ^ (2;). We 
treat the same three components as in the /esc = case: 
the truncated stellar spectrum, VJ^ , the free- free spectrum, 
(f, z, Zi), and the scattered Lya spectrum, l]^{v,z,Zi). 
The truncated stellar spectrum is given by eq. (051) , and the 
free- free spectrum is given by eq. (^ij). Similar to eq. (p6|), 
the Lya spectrum is 

l]/'°'{z,Zi) — qhya hv-Lya 0(j'Lya — V,z), (37) 



(38) 



where 

qi^ya{z, Zi) — (?Lya + 9LyQ • 

The total spectrum emitted by a Pop III star and the 
surrounding IGM from the formation of the star at Zi to 
the redshift at which the Universe is reionized, z = Zrcion, 
expressed in the rest frame of the star when it forms and 
normalized to the lifetime of the star, is 

ll{v,Zi) = 



1 

+ - 

T 





dt 


I dz 






'dz 



with 



[l^{u', 2, Zi) + z, Zi)] , (39) 



(40) 
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4.1 Model results 



A [/im] 

Figure 6. The cosmic infrared background from Pop III stars. 
The ordinate is observed frequency multipHed by observed spe- 
cific intensity. Shown are curves for star formation in haloes with 
virial temperatures above critical temperatures T^rit = 400 K 
(molecular-hydrogen cooling, solid lines) and 10^ K (atomic- 
hydrogen cooling, dashed lines), (a) /osc = 0. (b) /osc = 1. For 
both /esc cases, there are two sets of two curves: the upper set 
is for the ongoing star-formation model; the lower set is for the 
single-burst model. The properties of the curves are explained 
in §4.1; this figure shows z^nd = 7. The points show the excess 
CIRB, with 2(7 errors, and are described in §4.2. 



Because this spectrum is in the frame of the star when it 
forms, IGM emission at later times (and thus lower values 
of z) appears in the spectrum at higher frequencies than it 
is emitted at; this is because we blueshift all of the IGM 
radiation back to redshift Zi. The spectra are normalized 
by integrating the emission from Zi to z — 7 and dividing 
by the lifetime of the star, r. The resulting spectrum isn't 
meaningful for the observation of a single Pop III stellar 
source plus ionized IGM region, because the photons are 
emitted at different times, but it is a useful spectrum for 
computing the CIRB. The spectrum of a single source at 
redshift z < Zi (such that the star no longer contributes 
directly to the spectrum) is 



ll{v,z, 



(41) 



this spectrum is in the frame of the emission at redshift z. 

Figure |shows IT, l^{zi), Zj;^"(2i) and ll{zi) for Zi = 15. 
The sharp emission peak at 1216 A results from collisionally- 
excited Lya emission shortly after the star forms. The broad 
peak centered near 1000 A is Lya emission from recombina- 
tions. At 608 A the Lya spectrum is cut off as a consequence 
of integrating the IGM spectrum to 2: = 7; the cut-off loca- 
tion is 1216 A (H- «reion)/(l + Zi) . 



We now convolve the processed spectra of Section ^ with the 
star- formation histories of Section |^ to calculate the Pop III 
contribution to the CIRB. The background is evaluated with 



poo 


dt 


/ dz 




'0 


dl 



'^iz)jtiz) 
1 + z 



(42) 



(e.g., Peebles 1993). Here fobs is the observed frequency, 7^ 
is the observed specific intensity, iy{z) = (l+z)!^obs and j'^{z) 
is the comoving specific emission coefficient. Both the star- 
formation history and the assumed recombination history 
determine j'i^{z): 



(43) 



where li, is either I'l or ll, r = 2x 1 0^ yr is the fiducial m ain- 
sequence lifetime of a Pop III star (Bromm et al. 2001b) and 
tp{z) is either i/'on(z) or 'i/'burst(2), the SFR per comoving 
volume. 

Results for /esc = and 1 are shown in Fig. ^. There 
are a total of eight curves in the figure, generated by vary- 
ing each of three parameters over two values: /esc = or 
1; Tcrit = 400 or 10'' K; star-formation mode of ongoing 
or single-burst. Models with /esc ~ always produce more 
CIRB at all wavelengths than /esc ~ 1 models. Models with 
Terit = 400 K produce more CIRB than Tcrit = 10'' K mod- 
els at all wavelengths except near 1 fim for single-burst star- 
formation models. The ongoing star-formation models pro- 
duce more CIRB at all wavelengths than single-burst star- 
formation models, except near 1 pm for the /esc ~ case. 

The sharp edge at 1 /xm in all of the curves is a result of 
our sharp truncation of Pop III star formation at z^nd ~ 7- 
The edge occurs at the redshifted wavelength of Lya from 
stars at Zend, i.e., (1 -I- Zend)1216 A. The reason that the 
CIRB curves for Tcrit = 400 and 10* K have similar values 
at 1 fim for single-burst star-formation models is that the 
corresponding SFRs at 2 = 7 are similar (see Fig. ^). 

The "bump" in the CIRB curves from 0.7 fim to 1 
results from stellar emission between rest-frame 912 and 
1216 A (see Figs. | and ^) by sources at 2 = 7. From 
A = 1 /xm to roughly 2 /xm the Pop III CIRB is dominated 
by Lya emission from sources at 1 + z = A/1216 A. The in- 
fiection points of the curves, particularly clear in the /esc ~ 
panel, occur at the wavelength where the Pop III CIRB tran- 
sitions from being dominated by Lya emission to being dom- 
inated by continuum emission from sources at 2 = 2end- The 
continuum radiation for the /esc = spectrum is dominated 
by free- free emission (see Fig. ^ , whereas stellar continuum 
is important as well for the /esc ~ 1 continuum spectrum 
(see Fig. Consequently, the Pop III mid-IR background 
is determined predominantly by the choice of 2end rather 
than the SFR at extremely high redshifts. The Pop III con- 
tribution to the optical background (for 2end > 7) is due 
to the stellar "bump" described above and the high-energy 
tail of free-free emission, which is only significant for the 
/esc ~ spectrum; this is a consequence of the temperature 
of the gas, which is effectively much higher in the /esc = 
case than for the /esc ~ 1 case (see §3). 



4 COSMIC INFRARED BACKGROUND 
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4.2 Observational data 

The points with error bars in Fig. ^ show the difference 
between the total extragalactic background and the extra- 
galactic background due to resolved sources. At wavelengths, 
e.g., 0.3 /im, where the lower limit is an arrow, the con- 
tribution from resolved sources is sufficient to explain all 
of the background, at the 2cr level. At wavelengths where 
data points have lower limits plotted, e.g., 2.2 /im, the light 
from resolved sources is not enough to account for the mea- 
sured background. I t is t his unexplained excess, measured 
by Camb resy et al. ( 2001 ) at 1.25 wn and by Cambresy et 
al. ( |200l| ) and Wright & Johnson ( |200l|) at 2.2 pm, that we 
are trying to fit with our models, while conforming to upper 
limits at other wavelenths. 

Near-IR data points The triangles at 1.25 and 2.2 /xm 

(J and K band, respectively) correspo nd to the unexplained 
CIRB as measured by Cambresy et al. (2001): they measured 
the total extragalactic background with COB-E/DIRBE 
data, and subtracted from that the contribution due to re- 
solved galaxies measured by near-IR surveys. The squares 
at J and K (offset slightly in wavelength for clarity) 
are our computation of the unexplained CI RB fr om the 
CIRB measurement of Wright & Johnson (2001), using 
their measurement of the total extragalactic background 
from COBiS/DIRBE data (but using a different foreground 
model from Cambresy et al. 2001), and the same galaxy con- 
tribution subtraction and error pro pagat ion as Cambresy et 
al. ( |200l| ). The Wright & Johnson ( 200l| ) point at J band is 
consistent at the 2g level with no unexplained extragalac- 
tic background. All error bars are 2a. The 2a upper limit 
at 3.5 ^m (L band) represents a measurement of the to- 
tal ext ragalact ic background at that wavelength (Wright & 
Johnso:i 2001). SIRTF is expected to accurately measure 
the contribution of galaxies to the L-band CIRB, and thus 
allow a determination of the excess L-band CIRB (W. T. 
Reach, private communication). 

Other data points The upper limits at 0.3, 0.6, and 
0.8 are 2a upper limits on the unexplained optical extra- 
galactic background light (Bernstein, Freedman, & Madore 
2002). The open circle data points show the estimated un- 
explained optical extragalactic background, which is consis- 
tent with zero {2a) at all three bandpasses. Bernstein et 
al. (2002) measured the total optical extragalactic back- 
ground from HST images (utilising simultaneous ground- 
based spectra for absolute zodiacal-light subtraction cali- 
bration). From published number counts and their own care- 
ful photometry of the Hubble Deep Field observations, they 
find no evidence for convergence in the integrated light from 
galaxies; this suggests that deeper galaxy photometry will 
lower the upper limits. 

The upper limit at 15 jxm is our computation of 
the 2a upper limit on the unexplained mid-IR back- 
ground light, using the CIRB upper limits obtained by Re- 
nault et al. (2001) from gamma-ray observations, and the 
/SO/ISOCAM 15 /im galaxy counts reported by Elbaz et 
al. (2002). The open hexagonal data point shows the esti- 
mated excess CIRB, which is consistent with zero (2a). 



5 DISCUSSION 

Our models identify a very narrow range of parameter space 
in which Pop III stars may explain the excess CIRB: 

(i) The Pop III stars must be very massive, so that most 
of their energy is radiated in photons energetic enough to 
ionize hydrogen. 

(ii) Cooling in low-mass haloes is possible due to H2. 

(iii) Approximately 40 per cent of the eligible baryons in 
the star forming halo must be converted into Pop III stars. 

(iv) The escape fraction of ionizing photons from the neb- 
ula surrounding a Pop III star must be near zero, so that 
most of the ionizing photons are converted into Lya pho- 
tons. 

(v) Pop III star formation must begin by z ~ 25 and 
persist until 2;end — 7, so that the Lya emission (Fig. ^ 
extends through the observed J and K bands. 

(vi) Negative feedback effects must not inhibit Pop III 
star formation. In particular, we require that Pop III stars 
do not radiatively or mechanically destroy star-forming ma- 
terial, and that they they do not enrich their surroundings 
with a sufficient amount of metals to end Pop III star for- 
mation. 

The amplitudes of the curves in Fig. ^ scale simply 
with rj, and they do not change dramatically with small 
changes to the cosmological parameters. The characteristic 
breaks in the spectra are located at 912 A(l + Zcnd) and 
1216 A(l -I- Zond), where Zcnd is the low-redshift limit to 
Pop III star formation; if z^nd were much higher than 7, 
the Pop III CIRB peak would lie longward of the J band. 

If all of the above conditions are met, then all of the 
observed near-IR CIRB deficit is due to Pop III star for- 
mation. Future observations may demonstrate that other 
sources contribute significantly to the near-IR background; 
e.g., part of the deficit could be due to the faint wings of 
galaxies that are unaccounted for in current surveys (Totani 
et al. 2001). The parameter space explored in this paper 
easily accomodates lower values of the unexplained CIRB, 
through decreasing the efficiency, r). The shape of the Pop III 
CIRB spectrum in the near-IR is primarily determined by 
the shape of the star- formation rate as a function of redshift. 

The formation of very massive stars probably requires 
that the star-forming gas is not enriched with heavy ele- 
ments to a mass fraction Zcmt <^ 10^"^ Zq, where Zq is the 
solar value (Bromm et al. 2001a). We require a large frac- 
tion of the mass processed through Pop III stars to end up 
in massive black holes, so as not to overproduce metals com- 
pared to the observed levels in the high redshift IGM and 
to avoid enriching the primordial gas too quickly to Zait- 
Recently, Schneider et al. (2002) have pointed out the pos- 
sible problem that if all Pop III stars collapsed into massive 
black holes, the IGM would always have zero metallicity. It 
seems natural to assume, however, that a small fraction of 
the Pop III stellar mass, Qz ~ ZQb ~ epifiiii is ejected 
into the IGM through pair-instability supernovae, even if 
the majority of the mass is permanently locked up in black 
holes. Assuming a ratio flm/QB ~ 0.1, the required pair- 
instability fraction to produce the critical level of metallicity 
Zcrit would be Epi ~ 10~^ (see also Oh et al. 2001). 

We thus conclude that if Pop III stars explain a large 
fraction of the near-IR background, then almost all Pop III 
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mat eri al must end up in massive black holes. In that case we 



Agol E.. Kamio nkowski M., 2001, MNRAS, submitted 



would ] )redict a significant contribution, ^ 10 per cent, by 



pli/0109539) 



massive black holes to the total baryonic mass budget in the 
Universe. Although such a prediction may seem somewhat 
extreme, it will ultimately be tested by observations (e.g., 
Agol & Kamionkowski 2001; Schneider et al. 2002). We em- 
phasize that at present there are no observations that rule 
out such a scenario (e.g., Carr 1994). Moreover, there are 
good physical reasons to seriously consider the existence of 
massive black holes resulting from Pop III star formation 
(e.g., Bromm et al. 2002; Madau & Rees 2001; Schneider et 
al. 2002). 



In §2.4 we discussed the cumulative matter den- 

sity processed through Pop III stars. In Fig. ^ the horizontal 
lines labelled Hi, He I and Hen correspond to the fim re- 
quired to produce 10 ionizing photons per particle of each 
species (Bromm et al. 2001b). For the /esc ~ 1 model, the 
intersection of those lines with the Qui star curves is ex- 
pected to be closely related to the redshift of reionization. 
For the /esc = model, though. Pop III ionizing photons 
are all absorbed in the haloes in which they are emitted, 
and thus Pop III stars make no contribution to the reion- 
ization of the Universe. We discussed above that if Pop HI 
stars enriched their immediate environment with a sufficient 
amount of metals. Pop III stars would no longer form there 
(producing a SFR similar to our single-burst model) . If met- 
als from Pop III stars efficiently mixed throughout the Uni- 
verse, the era of Pop III stars might therefore come to an 
end (see Bromm et al. 2001a). 

A single Pop III star produces 
the Lya line, assuming /< 
in a flux of ~ 10 







esc — 0. At Z 
2^ —1 —2 

erg s cm 



M 



© 



10^'^ erg s-i M, 
: 7 that would result 
For a discussion of 
t he si ze of the Lya emitting region, see Loeb & Rybicki 
( 1999 ). For the /esc = 1 case the luminosity and size of the 
emitting region of a Pop HI star depend in more detail on 
redshift and IGM parameters. 

Future theoretical work should improve our understand- 
ing of the efficiency of Pop HI star formation, as well as the 
physical conditions of the haloes in which Pop III stars form. 
We finally expect the launch of NGST to open an obser- 
vational window into the earliest epochs of star formation, 
possibly including Pop HI stars, less than a decade from 
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